Power law tail in the radial growth probability 
distribution for DLA 

Peter Ossadnik and Jysoo Lee 
Hochstleistungsrechenzentrum (HLRZ) 

Forschungszentrum Jiilich GmbH 
Postfach 1913, W-5170 Jiilich, Germany 

February 1, 2008 

Abstract 

Using both analytic and numerical methods, we study the radial growth 
probability distribution P(r, M) for large scale off lattice diffusion lim- 
ited aggregation (DLA) clusters. If the form of P(r, M) is a Gaussian, 
we show analytically that the width £(M) of the distribution can not 
scale as the radius of gyration Rq of the cluster. We generate about 
1750 clusters of masses M up to 500, 000 particles, and calculate the 
distribution by sending 10 6 further random walkers for each cluster. 
We give strong support that the calculated distribution has a power 
law tail in the interior (r ~ 0) of the cluster, and can be described by a 
scaling Ansatz P(r, M) oc ^- • g ( r ~^ ^ , where g(x) denotes some scal- 
ing function which is centered around zero and has a width of order 
unity. The exponent a is determined to be ~ 2, which is now substan- 
tially smaller than values measured earlier. We show, by including the 
power-law tail, that the width can scale as Rg, i£ a > Df — 1. 

1 Introduction 

The growth of DLA [0 clusters can be described by a set of growth proba- 
bilities {p^ . Each of them describes the probability that site i is touched by 
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the next incoming particle. The determination of the pi and their distribution 
has gained much interest 0, [| [|, §]. Numerical calculations in this field 
were done by solving exactly the Laplace equation on a given DLA cluster 
[§, 0. Unfortunately these calculations are limited by computer resources 
to cluster masses around 50,000 particles. An alternative quantity, which 
is more accessible by large scale simulations, is the integrated radial growth 
probability P(r, M) . It describes the probability that the next incoming 
particle touches the cluster of mass M at a distance r from the seed. This 
quantity has been measured by Plischke and Racz || who studied the first 
two moments of the distribution and fitted a Gaussian behavior 

P(r,M) ocexp^-^^. (1) 

Using 4000 clusters of masses up to 2,500 particles they obtain a power law 
behavior of the center r$ and the width £ of the Gaussian 

r oc R G oc M u , £ oc M v> (2) 

with exponents v = 1/Df ~ 0.585, v' pa 0.48 , where Df denotes the fractal 
dimension of the clusters and Rq is the radius of gyration. Later, using 
larger scale simulations, Meakin and Sander || showed that the exponent 
v' approaches v with increasing cluster mass and hereby raised the question 
about the real behavior of the width £ . Thus, currently three different 
possibilities are being considered. 

1. The width scales with the same exponent v as Rq '■ oc M v . 

2. The width scales with a smaller exponent v' : £(M) oc M v \ v' <v 

3. The width scales with v but has logarithmic corrections: £(M) oc 
M7(ln(M)) /3 



Especially case (3) has attracted some interest [10|, [TT|. For an exponent 
(3 = 1/2 it implies a multiscaling behavior of the mass density M{x) oc r D<yX ^~ l 
where D(x) is a non constant function of x = r / Rq ■ 

In the following we will show that case (1) combined with the assumption 
of a Gaussian distribution leads to an unphysical singularity in the mass 
distribution in the limit of infinite M . This result is a consequence of the 
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fact that the Gaussian does not drop fast enough at the center of the cluster, 
which suggests that v' = v cannot be true. 

But on the other hand, we find numerically that a Gaussian is not a good 
description for the growth probabilities near the seed of the cluster, since for 
small r the tail behaves like a power law . Furthermore one has to realize 
that P(r, M) drops to zero for large r simply because of the finite size of the 
cluster. Usually such an effect would be taken into account by considering a 
finite size cutoff function. 

Here, we study the growth probabilities P(r, M) . We will show that 
the numerical data is only consistent with a Gaussian behavior around the 
maximum but that this description indeed fails in the small r tail. We will 
show that another type of behavior P(r, M) oc ^■g(j-jP L '*) where g{x) is a 
scaling function is consistent with our data as well around the maximum as 
in the tail. When this power law term is included in the distribution the 
previously mentioned singularity disappears. 

Thus, the organization of the paper is as follows. In section 2 we analyze 
the consequences of the different possible behaviors of £(M) on the mass 
density. In section 3 we study numerically the behavior of the growth prob- 
abilities and in section 4 we show the consequences of the power law tail on 
the mass density. 



2 Analytical results 

Following [7], we define two lengths, the radius r (M) and the width £(M), as 
r (M) = (r),£ 2 (M) = (r 2 )-(r) 2 , and r = \r\. Here, (f(r)) = J f(r)P(r,M)dr. 
Let us assume that P(r, M) can be described by a Gaussian distribution, 

P{r,M) = 1 exp[-(r-ro(M)) 2 /2£ 2 (M)]. (3) 
/2tt£ 2 (M) 



Here, we check whether the above form of P(r, M) is consistent with the 
growth of a fractal. We start from the relation 

N(r, M) = j P(r, M')dM', (4) 

where N(r, M)dr is the number of sites within [r, r + dr] in clusters of mass 
< M. Equation (4) can be derived from the fact that any site that lies within 
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[r, r + dr], in a cluster of mass M, has to be attached at the previous stage 
of growth (1 < M' < M) [9]. 

Case (1): Let us assume that r (M) = AM V and f(M) = BM U , where 
A, B are constants and we neglected other corrections. Using Eq. (4), we get 

N ( r > M ) = vW A M W^ e M-(r - AM'-f/{2B^M'^)]dM' 

= \/2ttB 2 eX P(~^2) A A? 717 eX P(~2B 2 M' 2 " + B^M'^dM 1 '. (5) 

Changing the variable of integration to x = M' ~ u Eq. (5) becomes 

The total mass contained in the [r, r + dr] shell is, by definition, N(r, M)dr. 
For a fixed value of r, iV(r, M) becomes larger for larger values of M. Es- 
pecially, in the M — > oo limit, the integral diverges due to the singularity at 
x = since l/i/ is larger than unity. This divergence is inconsistent with the 
well-established fractal description of DLA. 

Case (2): Consider the case £ = BM yv with < y < 1. The integral that 
corresponds to Eq. (5) becomes 

1 r M 1 

N(r,M) = -^==j i —eM-(r-AM>n 2 /(2B*M>*yn]dM>. (7) 
Changing the integration variable to x = M'~ v ', we get 

M = ^£-^ , " 1 " 1/ " exp (-^ 2( "" 1) -^ 2,+ ^" 1 )^ 

(8) 

In the M — > oo limit, Eq. (8) becomes 

N ( r > M ) = ^vW ■/? x^-^eM-^irx - Af]dx 

oc r 1 /"" 1 . (9) 

Since Df = l/u, Eq. (9) becomes N(r,M) ~ r D/_1 , consistent with the fact 
that the cluster is fractal with dimension Df. 
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Case (3): Finally, consider the case 



N(r, M) 



1 [M v^TM 7 



M 2B 



y { rM-»-Af dM , 



(10) 



V2^B 2 J i M< 



Changing the variable to x = M' u , Eq. (10) becomes 



N(r, M) 



1 



/ y/— lnx x 

m- v 



-l/u+(rx-A) 2 /(2uB 2 



] dx. 



(11) 



y/lixv^B 2 



Since the logarithmic correction to the width appears not only in the normal- 
ization of the Gaussian but also in its exponent, the behavior of the integrand 
is now substantially changed as compared to case (1). There, one observes 
a behavior x~ x l v ■ exp(— ax 2 + bx + c) (Eq. (6)), whilst here one obtains 
x -i/v+ax +bx+c_ rj^g integral in Eq. (11) can not be evaluated in a closed 
form, but we can still extract some of the properties. First, consider the 
M — > oo limit. The logarithmic term in the integrand is divergent at x = if 
the exponent of the power law term (A 2 — 2B 2 )/(2uB 2 ) is smaller than — 1. 
Therefore the whole integral is convergent only if A 2 > 2(1 — v)B 2 . If the 
integral in Eq. (11) is convergent, the convergent value can be estimated by 
treating the lnx term as a constant. The resulting form is consistent with 
the idea of multiscaling. 

The analysis presented above shows that if P(r, M) is Gaussian, the width 
of P(r, M), £(M), can not scale as Rq- The divergence in case (1) suggests 
that there is not enough "screening" at the inside of the cluster. Therefore, 
£(M) should increase slower than Rq as M is increased. 

In the next section, we numerically check whether the P(r, M) can well 
be described by a Gaussian. 

3 Numerical simulations 

To measure the growth probabilities P(x, M) numerically (here we scale 
x = r/R G ) we generate clusters of masses 10,000, 20,000, 50,000, 100,000, 
200,000 and 500,000 particles. For each mass we generate at least 150 dif- 
ferent clusters. Only for the 100,000 particle clusters we grow 1000 samples 
to have very good statistics and to see the details of the distribution. The 
growth probabilities are obtained in a static measurement: First we grow 
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each cluster to the full size and afterwards we probe its surface using an- 
other 10 6 random walkers. Using a fast algorithm [pj] , the generation of one 
100,000 particle cluster and the measurement of P(x, M) can be done within 
90 minutes on a Sun SPARC station 2. To obtain the P(x, M) we bin it for 
x G [0, 2.5] at 64 equidistant points and count the number of particles touch- 
ing the cluster at a radius between x and x + dx . Fig. [I] shows the growth 
probabilities for all cluster masses in linear and logarithmic scale. One has 
to notice the power law tail of P(x, M) at small x , which is a property that 
cannot be found for pure Gaussians. To smoothen out the noise we show in 
fig. the integrated distributions IP(x,M) = Jq P(x' , M) dx' . For small 
cluster masses up to 100,000 particles one determines a systematic decrease 
of the width. For the largest cluster masses no such a clear systematic be- 
havior can be seen. The curves for the 200,000 and 500,000 particle clusters 
lie on top of each other. The width £ can be studied quantitatively with two 
methods. 

The first one determines the width by measuring the first and second 



moments of the data £/Rg = y (x 2 ) — (x) 2 and the second consists in fitting 
a Gaussian to the data using a nonlinear fitting routine. Both methods give 
comparable results and one obtains 

^Moments x K G' M n9 \ 

£ rsr E> fl//— 0-050±0.003 

^Gauss K G ' 

We also check whether our data is consistent with case (3) and we find the 
behavior 

foments « /[n(107 . M) - (13) 



In fig. |3] we show the data of the width ^Moments/^ as a f unc tion of the 
cluster mass in semi logarithmic scale. The solid lines in this plot are the 
results obtained by fitting a power law and a 1/vhiM law to our data. For 
cluster masses up to 200,000 particles one seems to have good agreement 
with the assumption of a logarithmic correction. Thus, our data seems to be 
inconsistent with case (1) v' — v . For clusters up to 200,000 particles our 
data are consistent with cases (2) and (3). 

Here, we want to study in more detail the form of the inner part of 
P(x, M) . Since fig. [l| shows a clear power law behavior for small x , we 
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assume a simple scaling law 



x 



a 



X — Xq 



P(x,M) oc 



:f 



C/Rg 



(14) 



where f(x) is some scaling function. This type of behavior is used to stress 
the power law behavior of the small x tail. 

If one uses for f(x) too simple a function like a Fermi function which 
becomes constant for small x , one does not describe the full mass dependence 
of P(x, M) correctly: one measures only apparent large exponents a , which 
increase logarithmically with M (fig. |j). 

To take this mass dependence correctly into account we perform a data 
collapse to the scaling form (|i~4"|) . Here we use for the width £ / Rq and center 
Xq the values calculated from the moments of P(x, M) and vary the exponent 
a. The best collapse of our data, which is shown in fig. ([5]), we obtain for 
a value a = 2.0±0.4 . This value of a is much smaller than those obtained 
before and in 0. Thus, the previously measured large exponents are only 
"apparent" values. They are combinations of the power law tail and the mass 
dependence of the width of the scaling function. If one takes into account the 
mass dependence of the scaling function the a decrease very much. Anyhow, 
we still find a strong power law screening of the inner parts of the cluster. 

4 Consequences 

Having established the existence of the power-law tail of P(r, M) , we study 
how this modifies the results in Sect. |2|. Here we assume that P(r, M) can 
be described by eq. ( |HD with a normalization factor iVo . Let the average of 
f(x) be zero and the width be some constant. We first consider the case (1). 
Since P(r, M) has to be normalized, it is 




(15) 



Changing the integration variable to x = r/£ eq. fli~5|) becomes 




(16) 
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where xo = r /£. Since the integral is a constant, N Q should scale as {Rq/CT ■ 
The equation that corresponds to (5) is 



The part that causes a divergence in (5) is the M~( 1+a V D f term, since other 
terms in the integrand become constants for a fixed value of r . Therefore, 
the integral (|17|) diverges if the exponent (1 + a)/Df is less than unity, which 
is equivalent to a<Df — 1 . Since the numerical value of a is obviously larger 
than a = Df — 1 ~ 0.71 , our Ansatz is not anymore excluded by the self- 
consistency argument. One can show, in a similar way, that the existence of 
a power-law tail is consistent with the other cases (2) and (3). 



5 Conclusion 

We have studied the growth probabilities for large off lattice DLA clusters. 
We have shown that the usual assumption of a Gaussian behavior leads to 
unphysical singularities under the assumption that the width of the Gaussian 
scales with the same exponent as its center. We have shown that extensive 
numerical calculations seem to be inconsistent with the assumption v' = v . 
It seems that cases (2) or (3) have to be favored. Moreover, we showed that 
P(r, M) has a power law tail which is inconsistent with the Gaussian behavior 
that is usually assumed. Therefore we suggested another type of behavior 
P(r, M) oc ^/(^y^) which removes the unphysical singularity described 
above. This type of behavior describes the maximum of the measured growth 
probabilities as well as the tail. We showed that one explicitly has to take 
into account the dependence of the width of the distribution on the cluster 
mass in order to find the correct exponent a ~ 2. If one uses a simple ad 
hoc assumption do determine the inner tail of the distribution one mixes the 
"real" power law term and contributions from the scaling function. 

Unfortunately our data does not enable us to draw any conclusions about 
the behavior of the minimum growth probabilities, as studied for example 
in H]. The data we measured is only an average of the growth probabilities 
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within one "shell" of the cluster. Thus, we only have data about the spatial 
distribution of growth probabilities and not about the distribution n(P) of 
the growth probabilities themselves. A measurement of n(P) is for small P 
probably not feasible with a random walker method. 

We gratefully acknowledge stimulating discussions with H. J. Herrmann 
and D. Stauffer. 
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Figure 1: Growth probability distributions for different cluster masses in 
linear and double logarithmic scale. 



Figure 2: Integrated probability distributions. 



Figure 3: Scaling of the width of the distribution as obtained by measuring 
the moments of the distribution. The two solid lines denote the fits of a 
power law and a 1/lnM law to our data. 
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Figure 4: Behavior of the "apparent" exponents of the power law with varying 
cluster mass. This demonstrates, that the mass dependence is important in 
order to obtain the correct exponent. 



Figure 5: Plot of the scaling function g(y) oc P(£ • y + x , M) • £ • (£ • y + x ) a . 
The maximum of the distribution collapses as well as the power law tail. 
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